Electric-field induced polarization paths in Pb(Zri x Ti x )03 alloys 
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Properties of Pb(Zri_ a; Ti a :)03 (PZT) for compositions x near the morphotropic phase boundary and 
under an electric field are simulated using an ab-initio based approach. Applying an electric field 
of [111] orientation to tetragonal PZT (e.g., £=0.50) leads to the expected sequence of tetragonal, 
A-type monoclinic, and rhombohedral structures. However, the application of a field of orientation 
[001] to rhombohedral PZT (e.g., X—0A7) does not simply reverse this sequence. Instead, the system 
follows a complicated path involving also triclinic and C-type monoclinic structures. These latter 
phases are found to exhibit huge shear piezoelectric coefficients. 
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The existence of a morphotropic phase boundary 
(MPB) is an important feature of the phase dia- 
gram of many technologically interesting perovskite al- 
loys, such as [Pb(Zn 1/3 Nb 2/ 3)0 3 ](i_ :c )[PbTi03]. (PZN- 
PT), pb(Sci/ 2 Nbi /3 )0 3 k_ B )[PbTi()3]x (PSN-PT), and 
Pb(Zri_ K Ti x )03 (PZT) ft). In these materials, the MPB 
is the boundary in the temperature-composition plane 
separating the rhombohedral phase, in which the electric 
polarization is along the [111] direction, from the tetrag- 
onal structure, in which the polarization lies along the 
[001] direction. Applying an electric field to these alloys 
near the MPB generates an anomalously large change of 
strain. Very recent theoretical [0] and experimental jj| 
studies strongly suggest that such large piezoelectric re- 
sponse is driven by the electric-field induced rotation of 
the electrical polarization P, rather than an electric-field 
induced change of magnitude of this polarization. 

Despite these recent advances, many features related to 
the effects of electric fields on structural and piezoelectric 
properties of perovskite alloys near the MPB remain puz- 
zling. For instance, based on calculations performed on 
the simple BaTiG-3 system, Rcf. p proposed that apply- 
ing an electric-field along the pseudo-cubic [001] direction 
in rhombohedral PZN-PT induces a rotation of the po- 
larization occurring along the rhombohedral-tetragonal 
path (i.e., P continuously rotates from the [111] to [001] 
direction). On the other hand, Ref. J3| suggests a more 
complex mechanism in the same PZN-PT system in 
which the polarization first follows the rhombohedral- 
tetragonal path for small electric field, and then jumps to 
a new path connecting an orthorhombic to the tetragonal 
structure. However, due to experimental restrictions, the 
authors of Ref. || were only able to observe the second 
(and new) polarization path. As a result, the following 
questions are currently unanswered: If there is indeed 
a new polarization path at higher electric-field, what is 
the transitional mechanism allowing such a change of po- 



larization path? Is this transitional mechanism sudden 
or does it happen over a broad range of electric field? 
Another aspect that is poorly understood is the striking 
difference between the electric-field behavior on opposite 
sides of the MPB: applying an electric field to a tetrag- 
onal alloy leads to a much smaller change of strain than 
is observed for a rhombohedral composition jjjj. Fur- 
thermore, most studies investigating piezoelectricity in 
perovskite materials focused on the CZ33 piezoelectric co- 
efficient. As a result, other coefficients that may exhibit 
even larger values could have been overlooked. On may 
also wonder what is the consequence of the polarization 
paths proposed in Ref. J3| on piezoelectricity. 

In this Letter we address the problems summarized 
above by investigating the effect of an electric field on 
the structural and piezoelectric properties of PZT alloys 
near the MPB. Starting with a composition that is tetrag- 
onal at zero field, we find that applying an electric field 
along the [111] pseudo-cubic direction leads to a polar- 
ization that simply rotates in the (110) plane from [001] 
to [111], i.e., from the tetragonal (T) to the rhombohe- 
dral (R) structure, via an intermediate A-type monoclinic 
(Ma) structure Q. The resulting d 33 piezoelectric coeffi- 
cient is only large very close to the Ma~R transition. On 
the other hand, starting with a rhombohedral PZT com- 
position, we find that applying an electric field along the 
[001] pseudo-cubic direction does not simply generate the 
reverse transition sequence R-Ma~T. Instead, the polar- 
ization begins to follow this (TlO)-plane R-Ma~T path 
at low fields, but before reaching the T phase, the po- 
larization crosses over to join the (100)-plane O-Mq-T 
path connecting the orthorhombic (O) to the tetragonal 
structure via an intermediate C-type monoclinic (Mq) 
structure (as hinted for PZN-PT in Ref. [iY. (The no- 
tation for monoclinic phases is that of Ref. H , with Ma 
and Mc phases having polarization lying in the (110) and 
(100) planes, respectively.) The crossing from the Ma 
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to the Mc structure occurs via an intermediate triclinic 
(Tri) structure that is stable over a substantial electric- 
field range. The overall transition sequence is thus R- 
M A -Tri-M C ~T. We find that both the triclinic and M c 
phases exhibit huge shear d\5 piezoelectric coefficients. 

We use the numerical scheme proposed in Ref. Q, 
which involves the construction of an alloy effective 
Hamiltonian (i? e if) @ from first-principles calculations 
[p|-|lO[ and its use in Monte-Carlo simulations to compute 
finite-temperature properties of PZT alloys. The effect 
of an external static electric field E on physical prop- 
erties is included by adding a — P • E term in H e g (see 
Ref. pi). In the present study, the temperature is kept 
fixed at 50 K, and well-converged results are obtained 
with 10 x 10 x 10 simulation boxes (5000 atoms) |12|. We 
use atomic configurations mimicking maximal composi- 
tional disorder, consistent with experimental conditions 
p3| . Data from the simulations on the local mode vec- 
tors u (directly proportional to the electric polarization) 
and the homogeneous strain tensor are used within the 
correlation- function approach of Refs [p 14 to obtain 
the piezoelectric response. Up to 2 x10 s Monte-Carlo 
sweeps are first performed to equilibrate the system, and 
then 2 x 10 5 sweeps are used to get the various statistical 
averages. 

Previous studies ||[l2|,[l5) have demonstrated that our 
effective Hamiltonian approach gives a very good account 
of experimental and direct first-principles findings in per- 
ovskite alloys. In particular, it successfully predicts the 
existence of three low-temperature ferroelectric phases 
for Pb(Zri_ 2; Ti 2; )03 solid solutions near its MPB @]: a 
tetragonal (T) phase for larger x compositions (x > 49% 
at 50 K), a rhombohedral (R) phase for smaller x compo- 
sitions (x < 47.5% at 50 K), and the recently discovered 
monoclinic (Ma) phase Q in between. The polariza- 
tion is parallel to the pseudo-cubic [001], [111], or [vvl] 
(Q < v < 1) direction for the T, R, or Ma phase, respec- 
tively. The electrical polarization can thus be viewed as 
rotating in a (110) plane from [001] to [111] as the Ti 
composition decreases in the monoclinic phase (In 
the following, we adopt the convention that the x, y and 
z axes are chosen along the pseudo-cubic [100], [010] and 
[001] directions, respectively, and the piezoelectric tensor 
is represented in the orthonormal basis formed by ai = 
[100], a 2 = [010] and a 3 = [001].) 

We first consider a tetragonal PZT composition, 
a;=0.5, under an electric field applied along the [111] di- 
rection. The polarization of this alloy is initially along z 
when no electric field is applied. Figure 1(a) displays 
the cartesian components (u x , u y and u z ) of the su- 
percell average of the local mode vectors in this solid 
solution as a function of the magnitude of the electric 
field. It demonstrates that applying a small electric field 
along the [111] direction in Pb(Zr .5oTi . 50)03 perturbs 
the system in the expected way: the polarization rotates 
away from the [001] direction, as evidenced by the fact 
that u x and u v become non-zero. The strain tensor given 
by our simulations indicates that the resulting structure 
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FIG. 1. (a) Cartesian components u x , u y and u z of the 
local-mode vector, and (b) d$3 and dis piezoelectric coeffi- 
cients, as a function of the magnitude of the electric field 
applied along the pseudo-cubic [111] direction in disordered 
single crystals of Pb(Zro.5oTio. 50)03 at T=50K. Insert illus- 
trates the path followed by the polarization. 



is the monoclinic Ma structure. As the magnitude of 
the electric field along the [111] direction increases, u z 
decreases, while u x and u y increase and remain equal. 
The polarization in this Ma structure thus rotates to- 
wards the [111] direction within the (110) plane. For a 
sufficiently large electric field (~110 x 10 6 V/m), u z be- 
comes equal to u x and u y , si gna ling a phase transition 
to a rhombohedral structure |16|] with the polarization 
aligned along the [111] direction. The path followed by 
the polarization in going from the T via the Ma to the R 
structure with increasing field is illustrated in the insert 
of Fig. 1(a). 

Figure 1(b) shows the resulting behavior of the CZ33 
and di5 piezoelectric coefficients. We decided to focus 
on these two coefficients because we numerically found 
that they are the largest ones among all the 18 piezoelec- 
tric constants. One can note that di§ is larger than c?33 
in the T structure (and in most of the range of stabil- 
ity of the Ma structure). This is consistent with recent 
measurements revealing that the piezoelectric elongation 
of the tetragonal unit cell of PZT does not occur along 
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the polar [001] direction (Tt]]. The large value of di$ also 
explains the strong (averaged) piezoelectric response ob- 
served in ceramic tetragonal samples 0j. One can also 
see that ^33 is predicted to be larger in the R phase than 
in the T and Ma structures. This prediction is consistent 
with the findings of Refs. [PJl^| that the ^33 coefficient 
of rhombohedral materials can be very large along the 
[001] direction when the polar direction is oriented along 
the pseudo-cubic [111] direction. Figure 1(b) also re- 
veals that (I33 adopts values larger than 1,000 pC/N in a 
narrow electric-field range centered around the MA-to-R 
transition. Note that G?i5 is also enhanced in this region, 
but has much smaller values than ^33. 

We now investigate the effect of an electric field applied 
along the [001] direction to a rhombohedral PZT com- 
position, Pb(Zro.53Tio. 47)03. The polarization of this 
alloy starts along the pseudo-cubic [111] direction (i.e, 
) when no external electric field is present. 
Figure 2(a) reveals an unusual behavior of the polariza- 
tion and a resulting rich phase diagram. As expected, 
applying a small electric field along the [001] direction of 
Pb(Zr .53Ti . 47)03 perturbs the R structure to the Ma 
structure, for which < u x = u y < u z . As the magnitude 
of the electric field increases, the polarization rotates to- 
wards the [001] direction, with u z increasing while u x 
and u y are decreasing. However, as shown in Fig. 2(a), 
Pb(Zi'o.53Tio. 47)03 then adopts two behaviors that do 
not occur in the Pb(Zro.5oTio. 50)03 alloy. First, for an 
electric field ranging between 50 and 88 x 10 6 V/m, u x , 
u y and u z are all non-zero and different from each other, 
and the strain tensor identifies this as a triclinic (Tri) 
phase. Second, for an electric field ranging between 88 
and 110 x 10 6 V/m, u x vanishes while u v and u z remain 
zero and different from each other. This characterizes a 



second type of monoclinic phase denoted Mc [5|, |l9| . pi0[ . 
At larger electric field, Pb(Zr . 53^0.47)03 finally trans- 
forms from this monoclinic (Mc) phase into the tetrag- 
onal (T) phase, for which the polarization is along the 
[001] direction [u z is then the only nonzero component). 

The insert of Fig. 2(b) summarizes the path followed 
by the polarization in Pb(Zro. 53^0.47)03: P first rotates 
in the (110) plane from [111] to [vvl] (R to M A ), then 
switches via [wvl] (Tri) to the (100) plane, where it fin- 
ishes by rotating from [Ovl] (M c ) to [001] (T). That is, 
the polarization starts along the R-Ma~T path and then 
switches to the O-Mc-T path to complete the transfor- 
mation to the T phase. 

The electric-field induced R-Ma~T transformation 
path has been observed exp erimentally in nominally 
rhombohedral PZT alloys WA , but the unavailability of 
single-crystal PZT complicates the experimental deter- 
mination of the polarization paths depicted in Fig. 2(a). 
Furthermore, we find that the change of path occurs at 
electric fields that are high compared to those that can 
be used experimentally (our calculations predict that an 
electric field of ~ 100 x 10 6 V/m is needed to see the 
transitions to the Tri and Mc phases, while the highest 
field applied experimentally in Ref. ft is 2 x 10 6 V/m). 
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FIG. 2. (a) Cartesian components u x , u y and u z of the 
local-mode vector, and (b) and dis piezoelectric coeffi- 
cients, as a function of the magnitude of the electric field 
applied along the pseudo-cubic [001] direction in disordered 
single crystals of Pb(Zro. 53^0.47)03 at T=50K. Insert illus- 
trates the path followed by the polarization. 



However, we find that the transformations displayed in 
Fig. 2(a) all occur at smaller electric fields in PZT alloys 
for which the Ti composition x is closer to the MPB: 
increasing the Ti composition in PZT from 46% to 47% 
leads to a considerable decrease of the critical electric 
field, giving rise to the Tri phase at ~ 50 instead of 
f50xl0 6 V/m, and to the Mc phase at ~ 88 instead 
of 237 x f0 6 V/m. Thus a practical observation of the 
change of polarization path depicted in Fig. 2(a) would 
probably require the use of a single-crystal PZT alloy 
with a composition lying as close as possible to the MPB. 

While the R-Ma~T path of the polarization has re- 
cently been predicted to occur in (zero-field rhombohe- 
dral) BaTi03 when applying an [001]-oriented electric 
field, no change of path has been predicted to occur in 
this material even at very high field ^jn}. Similarly, 
the O-Mc-T path has recently been observed in (zero- 
field rhombohedral) 92%-PbZn 1/3 Nb 2 /30 3 8%-PbTi0 3 
under an [001]-oriented electric field []3|. However, be- 
cause of some experimental limitations, the R-Ma~T 
path and the transitional triclinic phase, that probably 
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both occur at small fields, have yet to be observed in this 
material. To our knowledge, the present study is thus 
the first demonstrating the existence of these two differ- 
ent polarization paths, and the transitional phase con- 
necting them, in a cubic perovskite system under electric 
field @. 

Furthermore, Fig. 2(b) clearly shows that these un- 
usual polarization paths have a drastic effect on the g?33 
and c?i5 piezoelectric coefficients of Pb(Zr . 53Ti . 47)03. 
For instance, Fig. 2(b) reveals that the (I33 coefficient 
suddenly peaks at a small field, when Pb(Zr . 53^0.47)03 
adopts the Ma structure. This peak corresponds to a 
sudden increase of the strain 773. ^33 then decreases with 
a further increase of the magnitude of the electric field, 
implying that the strain 773 then adopts a flatter strain- 
vs. -field slope. Interestingly, both the sudden increase 
and the flatter behavior of the strain have been observed 
in rhombohedral PZN-PT 0| and were predicted to oc- 
cur in rhombohedral BaTi03 0] when applying an [001]- 
oriented electric field. 

The most striking result of Fig. 2(b) is the huge en- 
hancement of the shear d 15 piezoelectric coefficient as- 
sociated with the change of polarization path: d\§ in- 
creases from 1,000 to 8,000 pC/N with increasing electric 
field within the transitional Tri phase. Once the polar- 
ization is along the O-Mc-T path (i.e., when the system 
has transformed to the Mc phase), dis decreases with a 
further increase of the electric field However, dxs re- 
mains large, even in the tetragonal phase under very high 
electric field. As a result, tfx5 is larger than 2,000pC/N 
for a wide range of electric fields, e.g., ~ 50 x 10 6 V/m. 

In summary, we have used the first-principles derived 
computational scheme proposed in Ref. || to study 
the structural and piezoelectric properties of disordered 
Pb(Zri_ x Ti K )03 solid solutions near the MPB under an 
electric field. Our main results are: (1) the response of 
the polarization to an electric field is qualitatively dif- 
ferent in the tetragonal and rhombohedral phases; (2) 
specifically, rhombohedral alloys near the MPB exhibit 
an unusual switching of the polarization path under ap- 
plied field; (3) a large piezoelectric response can result 
from such a path switch; and (4) the shear piezoelec- 
tric coefficients (e.g., di$) could be exploited to achieve a 
drastic improvement in response in piezoelectric devices. 
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